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Abstract. The Mahonian statistic is the number of inversions in a 
permutation of a multiset with a; elements of type i, 1 < i < m. The 
counting function for this statistic is the q analog of the multinomial 
coefficient (™^ + + a am ), and the probability generating function is the 
normalization of the latter. We give two proofs that the distribution 
is asymptotically normal. The first is computer- assisted, based on the 
method of moments. The Maple package MahonianStat, available from 
the webpage of this article, can be used by the reader to perform exper- 
iments and calculations. Our second proof uses characteristic functions. 
We then take up the study of a local limit theorem to accompany our 
central limit theorem. Here our result is less general, and we must be 
content with a conjecture about further work. Our local limit theo- 
rem permits us to conclude that the coemecients of the g-multinomial 
are log-concave, provided one stays near the center (where the largest 
coefficients reside.) 



1. Introduction 

The most important discrete probability distribution, by far, is the Bino- 
mial distribution, B(n,p) for which we know everything explicitly, ¥(X = i) 
— p) n ~*), the probability generating function ((pt + (1 — p)) n ), the 
moment generating function {{pe l + 1 — p) n ), etc. etc. Most importantly, it 
is asymptotically normal, which means that the normalized random variable 

X n -np 

y/np(l -p) 

tends to the standard Normal distribution N(0, 1), as n — > oo. 

Another important discrete distribution function is the Mahonian distri- 
bution, defined on the set of permutations on n objects, and describing, 
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inter-alia, the random variable "number of inversions" . (Recall that an in- 
version in a permutation m,...,ir n is a pair 1 < i < j < n such that 
t^i > TTj)- Let us ca ^ this random variable M n . The probability generating 
function, due to Netto, is given explicitly by: 

The formula (1.1) has a simple probabilistic interpretation (see Feller's 
account in [3, Section X.6]): If Yj is the number of i with 1 < i < j and 
7Tj > nj, then 

M n = Y 1 + • • • + Y n , (1.2) 

and Yi,...,Y n are independent random variables and Yj is uniformly dis- 
tributed on {0, . . . ,j — 1}, as is easily seen by constructing ir by inserting 
1, . . . ,n in this order at random positions; thus Yj has probability gener- 
ating function (1 — q^)/(j{l — <?))• It follows from (1.1) or (1.2) by simple 
calculations that the Mahonian distribution has mean and variance 

EM n = ^l, (1.3) 

VarM = n ( n ~ 1 )( 2n + 5 ) = 2n 3 + 3n 2 -5n 

72 72 1 ' J 

Even though there is no explicit expression for the coefficients themselves 
(i.e. for the exact probabilitity that a permutation of n objects would have 
a certain number of inversions), it is a classical result (see [3, Section X.6]), 
that follows from an extended form of the Central Limit Theorem, that the 
normalized version 

M n - n(n - l)/4 

y/(2n 3 +3n 2 - 5n)/72 ' 

tends to N(0, 1), as n — > oo. So this sequence of probability distributions, 
too, is asymptotically normal. 

But what about words, also known as multi-set permutations?. Permu- 
tations on n objects can be viewed as words in the alphabet {1,2, . . . ,n}, 
where each letter shows up exactly once. But what if we allow repetitions? 
I.e., we consider all words with a\ occurrences of 1, a2 occurrences of 2, . . . , 
a m occurrences of m. (We assume throughout that m > 2 and each a,j > 1.) 
We all know that the number of such words is the multinomial coefficient 

a± + ■ ■ ■ + a r , 

&1 j • • • j Q"m 

and many of us also know that the number of such words with exactly k 
inversions is the coefficient of q k in the g-analog of the multinomial coefficient 

ai H h a m \ [ai H h a m ]\ . . 



ai, . . . ,a m J " [ai]! • • • 
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where [n]\ := [1][2] • • • [n], and [n] := (1 — q n )/(l — q); see [1, Theorem 3.6]. 
Assuming that all words are equally likely (the uniform distribution), the 
probability generating function is thus 

p (a) . TOW)-n^r +am (W) F ai+ ... +am {q) 

ai '-' amK9) ■ («! + ••• + a m )\ ut=i n-ii(i - <?) F *Ai) • • • ^(9) ' 

(1.6) 

Indeed, this can be seen as follows. Let M aiv .. )0m denote the number of 
inversions in a random word. If we distinguish the ctj occurrences of i by 
adding different fractional parts, in random order, the number of inversions 
will increase by Zj, say, with the same distribution as M ai ; further M 0lj ... j0m 
and Zi, . . . , Z m are independent. On the other hand, M aiv .. j<lm +Zi+- • -+Z m 
has the same distribution as M aiH |_ am . Hence, 

^,...^(9)^(9) ■■■^a m (9) = J F 01 +...+a m (?), (1-7) 

which is (1.6). 

By (1.6), we further have the factorization 

m 

F ai ,..,a m (q) = l[F Aj _ 1:aj ( q ), (1.8) 

J'=2 

where A,- := Gq + ■ ■ ■ + a^, which reduces the general case to the two- letter 
case. 

Note that (1.6) shows that the distribution of M 0lj .„ j<lm is invariant if we 
permute a±, . . . , a m ; a symmetry which is not obvious from the definition. 

Remark 1.1. The twodetter case is particularly interesting, since the un- 
normalized generating function 

+ b \ F / x _ (1 — q a+b ){l - q a+b ' 1 ) ■■■(!- q a+1 ) _ [0 + 6]! 

(the g-binomial coefficient in (1.5)) is the same as the generating function 
for the set of integer-partitions with largest part < a and < b parts, in other 
words the set of integer-partitions whose Ferrers diagram lies inside an a 
by b rectangle, where the random variable is the "number of dots" (i.e. the 
integer being partitioned). In other words, the number of such partitions of 
an integer n equal the number of words of a l's and b 2's with n inversions. 
See Andrews [1, Section 3.4]. 

It is easy to see that the mean of M 01j ... j0m is 

(x(ai, . . .,a m ) := EM olr . |(lm = e 2 (ai, . . . ,a m )/2 

(here efc(ai, . . . , a m ) is the degree k elementary symmetric function), so con- 
sidering the shifted random variable M 0lj ... j0m — /u(ai, . . . , a m ), "number of 
inversions minus the mean" , we get that the probability generating function 
is 

Q ( a ) ._ -/i(ai,...,a m ) p , \ _ Fai,...,am(Q) /-. q-, 

<^ai,...,a m W — q fa u ...,a m W ~ e2 ( aii ..., 0fc )/ 2 ^• y ' ) 
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By computing (q(qG)')' and plugging-in q = 1, or from (1.7) and (1.3)-(1.4), 
it is easy to see that the variance a 2 := Var M ai) ... )<lm is 

a 2 = (C1 + 1)62 " 63 . (1.10) 
12 v ; 

(By a we mean a(a±, . . . , a rn ) and we omit the arguments (ai, . . . , a m ) from 
the e^'s.) 

Let A 7 " := ei = ai + ■ ■ ■ + a m , the length of the random word, and let 
a* := maxj a,j and iV* := N — a*. 

One main result of the present article is: 

Theorem 1.2. Consider the random variable, M 0l) ... jtlm , "number of inver- 
sions", on the (uniform) sample space of words with a\ l's, ci2 . . . , a m 
m's. For any sequence of sequences (ai,...,a m ) = (a^ , . . . , ) suc/i 
i/iat iV* := A?" — a* — > oo, i/te sequence of normalized random variables 
v __ M au ..., am -/x(ai,... 

cr(ai, . . . ,a m J 

tends to the standard normal distribution JV(0, 1), as ^ — > oo. 

Theorem 1.2 includes both the case when m > 2 is fixed, and the case 
when m — > oo. If m is fixed and a\ > ai > • • • > a m , as may be assumed 
by symmetry, then the condition TV* — > oo is equivalent to a2 — > oo. In the 
case m — > oo, the assumption iV* — > oo is redundant, because N* > m — 1. 

Remark 1.3. The condition iV* — > oo is also necessary for asymptotic 
normality, see Section 5. 

We give a short proof of this result using characteristic functions in Sec- 
tion 3. We give first in Section 2 another proof (at least of a special case) 
that is computer- assisted, using the Maple package MahonianStat available 
from the webpage of this article: 

http : //www . math . rutgers . edu/~zeilberg/mamarim/mamarimhtml/mahon . html, 

where one can also find sample input and output. This first proof uses the 
method of moments. 

We conjecture that Theorem 1.2 can be refined to a local limit theorem 
as follows: 

Conjecture 1.4. Uniformly for all a±, . . . ,a m and all integers k, 

^ -» = fc ) = 7k(^" , ' ,3/,2 ° 2,+0 (^))- ('■") 



. 27TO- 

We have not been able to prove this conjecture in full generality, but we 
prove it under additional hypotheses on a\, . . . , a m in Section 4. 

2. A COMPUTER-INSPIRED PROOF 

We assume for simplicity that m is fixed, and that (a±, . . . , a m ) = (to,®, . . . , ta [ 
for some fixed a°, . . . , and t — > oo. 
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We discover and prove the leading term in the asymptotic expansion, 
in t, for an arbitrary 2r-th moment, for the normalized random variable 
X ai ,...,a m = (M 01j ... j0m — fJ,)/(T, and show that it converges to the moment 
fi 2r ='(2r)!/(2 r r!) of Af(0, 1), for every r. 

For the sake of exposition, we will only treat in detail the two-letter 
case, where we can find explicit expressions for the asymptotics of the 2r-th 
moment of M aii(l2 — fi, for a± = ta, a2 = tb with symbolic a, b, t and r to any 
desired (specific) order s (i.e. the leading coefficient t 3r as well as the terms 
involving t 3r_1 , . . . , t 3r ~ s ). A modified argument works for the general case, 
but we can only find the leading term, i.e. that 

a a r:=E(Jr ai> ... >am ) a - = M + 0(t -i) . 

Of course the odd moments are all zero, since the distribution of M 01j ... j0m 
is symmetric about fx. 

In the two-letter case, the mean of M a ^ is simply ab/2, so the probability 
generating function for M a ^ — fi is, see (1.6), 

_ F a , b (q) _ a\b\(l - q a + b )(l - q^ 1 ) •••(!- q^ 1 ) 
a ' b(q) q ab / 2 q ab / 2 (a + 6)!(1 - q b )(l - q^ 1 ) ■ ■ ■ (1 - q 1 ) ' 
Taking ratios, we have: 

G a ,b(Q) «(1 - Q a+b ) (2 ^ 

G a . hb (q) q b /i{a + b){l-q«) ' K ' ) 

Recall that the binomial moments B r := E,[( Ma ^~ f1 ^] are the Taylor coef- 
ficients of the probability generating function (in our case G a ,b(q)) around 
(7 = 1. Writing q = 1 + z, we have 

oo 

G a , b (l + z)=J2 B r(a,b)z r . 

r=0 

Note that Bq(cl, b) = 1 and B\(a,b) = 0. Let us call the expression on the 
right side of (2.1), with q replaced by 1 + z, P(a, b, z): 

l ' ' >' (l + Z )6/2(o + 6)(l-(l + z)«) ' 

Maple can easily expand P(a, b, z) to any desired power of z, It starts out 
with 

P(a, b, z) = 1 + ^ (2 a + b) bz 2 - ± (2 a + 6) bz 3 

— (8 a 3 - 8 a 2 b - 12 ab 2 - 3 6 3 - 440 a - 220 b) bz 4 + . . . 

5760 V ; 

note that the coefficients of all the powers of z are polynomials in (a, b). 
So let us write 



P(a,b,z) = ^2p i (a,b)z t , 



i=0 
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where Pi(a, b) are certain polynomials that Maple can compute for any i, no 
matter how big. 

Looking at the recurrence 

G a , b (l + z) = P(a,b,z)G a - 1>b (l + z) , 

and comparing coefficients of z r on both sides, we get 

r 

B r (a, b) -B r (a- 1,6) = ^B r _ s (a - l,b)p s {a,b) . (2.2) 

s=l 

Assuming that we already know the polynomials -B r _i(a, b), l? r _2(a, b), . . . , 
So (a, 6), the left side is a certain specific polynomial in a and b, that Maple 
can easily compute, and then B r (a,b) is simply the indefinite sum of that 
polynomial, that Maple can do just as easily. So (2.2) enables us to get 
explicit expressions for the binomial moments B r (a,b) for any (numeric) r. 

But what about the general (symbolic) r? It is too much to hope for the 
full expression, but we can easily conjecture as many leading terms as we 
wish. We first conjecture, and then immediately prove by induction, that 
for r > 1 

„ , „ 1 fab(a + b)\ r , 

r>2r(a, b) = — - I — I + lower order terms 

n , ,n - 1 fab{a + b)\ r , 

B 2r +i{a, b) = — - — I — j + lower order terms , 

where we can conjecture (by fitting polynomials in (a, b) to the data obtained 
from the numerical r's) any (finite, specific) number of terms. 

Once we have asymptotics, to any desired order, for the binomial mo- 
ments, we can easily compute the moments fj, r (a, b) of M a ^ — fx themselves, 
for any desired specific r and asymptotically, to any desired order. We do 
that by using the expressions of the powers as linear combination of falling- 
factorials (or equivalently binomials) in terms of Stirling numbers of the sec- 
ond kind, S(n,k). Note that for the asymptotic expressions to any desired 
order, we can still do it symbolically, since for any specific m, S(n,n — m) 
is a polynomial in n (that Maple can easily compute, symbolically, as a 
polynomial in n). In particular, the variance is: 

2 . ab(a + b+l) 
a = fJl2 (a,b) = — , 

in accordance with (1.10). In general we have /Lt2 r +i(a, b) = 0, of course, 
and the six leading terms of fi2r( a t, bt) can be found in the webpage of this 
article. From this, Maple finds that a2 r {at,bt) := p,2 r (at,bi)/ fj,2{at,bt) r are 
given asymptotically (for fixed a, b and t — > oo) by: 

(2r)\ ( r(r - 1) (b 2 + ab + a 2 ) l\ , „ 
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In particular, as t — ► oo, they converge to the famous moments of jV(0, 1). 
QED. 

2.1. The general case. To merely prove asymptotic normality, one does 
not need a computer, since we only need the leading terms. The above proof 
can be easily adapted to the general case (a±, . . . , a m ) = (to,®, . . . , ta^)- One 
simply uses induction on m, the number of different letters. 

2.2. The Maple package MahonianStat. The Maple package MahonianStat, 

accompanying this article, has lots of features, that the readers can ex- 
plore at their leisure. Once downloaded into a directory, one goes into a 
Maple session, and types read MahonianStat ; . To get a list of the main 
procedures, type: ezra();. To get help with a specific procedure, type 
ezra(ProcedureName) ;. Let us just mention some of the more important 
procedures. 

AsyAlphaW2tS(r,a,b,t,s): inputs symbols r,a,b,t and a positive inte- 
ger s, and outputs the asymptotic expansion, to order s, for ol^t (=A*2r/ 7*2 ) 

ithMomWktE(r,e,t): the r-th moment about the mean of the number 
of inversions of ait l's, . . . , a m t m's in terms of the elementary symmetric 
functions, in ai, . . . , a m . Here r is a specific (numeric) positive integer, but 
e and t are symbolic. 

AppxWk(L , x) : Using the asymptotics implied by the asymptotic normality 
of the (normalized) random variable under consideration, finds an approxi- 
mate value for the number of words with L[l] l's, L[2] 2's, . . . , L[m] m's with 
exactly x inversions. For example, try: AppxWk( [100,100,100] ,15000) ; 

For the two-lettered case, one can get better approximations, by procedure 
BetterAppxW2, that uses improved limit-distributions, using more terms in 
the probability density function. 

The webpage of this article has some sample input and output. 



We have an exact formula (1.10) for the variance a 2 of M ai) ... j0m . We first 
show that a 2 is always of the order Q(N 2 N*). 

Lemma 3.1. For any a\, . . . ,a m , 



Proof. For the upper bounds we assume, by symmetry, that a± > ■ ■ ■ > a. 
Then a* = a\ and 



3. A GENERAL PROOF OF THEOREM 1.2 



N 2 N* 

< 

36 ~ 



(N + 1)NN* N 2 N* 
12 ~ 6 




i=2 



3=3 



i=2 



Since e\ = N, (1.10) yields the upper bounds. 
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For the lower bound, we first observe that 2e2ei — 6es > (since this 
difference can be written as a sum of certain ajCikai). Hence ez < eie2/3 
and (1.10) yields 



12a 2 > eie 2 - e 3 > f eie 2 . 



Further, 



2e 2 = a j( N ~ a i) ^ a i( N ~ a *) = NN *> 

3=1 3=1 

and the lower bound follows. □ 
Proof of Theorem 1.2. From (1.6) follows the identity 

F - ft (^ + ^-l)/(i(n 1+ j)^) 



By Taylor's series 



\og^—^ = z/2 + z 2 /24 + 0(z A ), \z\<l, 

and we substitute this expansion into the identity (3.1) to conclude: 

F ni ,n 2 (e w ) = exp (inm 2 0/2 - n x n 2 (n x + n 2 + l)0 2 /2A + O(n 2 n 4 4 )) , 

(3.2) 

uniformly for ni > n 2 > 1 and \0\ < {n\ + n 2 ) _1 . 

We use the factorization (1.8). By symmetry, we may assume a\ > a 2 > 
■ ■ ■ > «m, and then > a^-i > aj for each j. Thus (3.2) yields, uniformly 
for q = e w with \0\ < N~\ 

m 

F ai ,...,a m (q) = Y[F Aj _ 1>aj (q) 

3=2 

(rn 
J^(iAj_iaj-0/2 - Aj^cijiAj + l)# 2 /24 + O( aj A 4 _i0 4 )) 

Here, the sums of the coefficients of 6 and 9 2 are easily evaluated, but we 
do not have to do that since they have to equal i/j, and —a 2 /2, respectively. 
Further, 

m m 

A)_ x a 3 <N 4 Y, *3 = N 4 N*. (3.3) 

3=2 j=2 

Consequently, if \9\ < N^ 1 , 

F au ...,a m (e w ) = exp(i/z0 - a 2 9 2 /2 + O(N 4 N*0 4 )) (3.4) 

and, by (1.9), 

G ai ,..., am (e w ) = exp(-a 2 2 /2 + 0(N 4 NJ 4 )). (3.5) 
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Let 6 = tja. For any fixed t, by Lemma 3.1, 

\Nt/a\ = 0(N; l/2 ) =o(1), 
so \9\ < iV -1 if v is large enough. Hence, by (3.5) and Lemma 3.1, 

G au ..., am (^) = ex p(-? + °(^f)) =ex P H 2 /2 + (l)), 
and Theorem 1.2 follows by the continuity theorem [4, Theorem XV. 3. 2]. □ 

4. The local limit theorem 

"If one can prove a central limit theorem for a sequence a n (k) of numbers 
arising in enumeration, then one has a qualitative feel for their behavior. 
A local limit theorem is better because it provides asymptotic information 
about a n (k) . . . ," [2]. In this section we prove that the relation (1.11) holds 
uniformly over certain very general, albeit not unrestricted, sets of tuples 
a = (ai, . . . , a m ). The exact statement is given below in Theorem 4.5. 

As explained in Bender [2], there are two standard conditions for passage 
from a central to a local limit theorem: (1) if the sequence in question is 
unimodal, then one has a local limit theorem for n in the set {|n — jj\ > ea}, 
e > 0; (2) if the sequence in question is log-concave, then one has a local 
limit theorem for all n. Our sequence, the coefficients of the ^-multinomial, 
is in fact unimodal, as first shown by Schur [6] using invariant theory, and 
later by O'Hara [5] using combinatorics. Unfortunately, the ensuing local 
limit theorem fails to cover the most interesting coefficients, the largest ones, 
near the mean [i. However, our polynomials are manifestly not log-concave 
as is seen by inspecting the first three coefficients (assuming n\,n2 > 2) 

( ni + n2 ) =l + , + 2, 2 + ... . 

V "1 J q 

The question arises might the coefficients be log-concave near the mean, and 
here is a small table of empirical values: (c[j] = [qi] ( 2 ™) q ) 

n (c[n 2 /2-l}) 2 -c[n 2 /2] x c[n 2 /2 - 2] 



2 -1 

4 -7 

6 -165 

8 -1529 

10 44160 

12 7715737 

14 905559058 

16 101507214165 

18 11955335854893 

20 1501943866215277 



Based on this scant evidence, we speculate that some sort of log-concavity 
theorem is true, but that its proper statement is complicated by describing 
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the appropriate range of a and j. Thus, we use neither of the two standard 
methods mentioned above for proving our local limit theorem. (Later, we 
shall see that our theorem has implications for log-concavity.) Instead, we 
use another standard method, direct integration (Fourier inversion) of the 
characteristic function, or equivalently of the probability generating function 
F(q) for q = e ld on the unit circle. We begin with one such estimate for 
rather small 0. 

Lemma 4.1. There exists a constant r > such that for any a\, . . . ,a m 

and \0\ < t/N, 

|i ? a 1 ,... > a m (e w )| = |G oll ... >am (e w )|<e-^/ 4 . 
Proof. Suppose that < \0\ < t/N. Then, using Lemma 3.1, 



N A N if A N 2 N*t 2 



^ <"^-<36r*, 

so if r is chosen small enough, the error term O{N A N*0 A ) in (3.4) and (3.5) 
is < a 2 2 /4, and thus the result follows from (3.5). □ 

We let in the sequel r denote this constant. We may assume < r < 1. 
Lemma 4.2. Uniformly, for all a\, . . . ,a m and all integers k, 



nM ai ,..., am = k)- 1 e -(^) 2 /(20 

V27TCT 



Proof. For any integer k, 

P(M ai ,..., am = fe)-^ e -(^) 2 /(^) 

V27T<T 



= i- f F au ..., am ^)e-- M d0 - i- f e-^e-^de 

= ±f {G ai ,.., am ^)-e^/Ae'^Ue 
2vr J\e\<r/N K ' 

+ I F ai _ am (e W )e- 

Z7T Jr/N<\e\<7T 

1 e-^^e-^-^dO 



~ike de 



2vr J\o\> T /N 
= : h+h + h- 

By (3.5) and the inequality \e w — 1| < \w\ max(l, \e w \) we find for \0\ < 
t/N, using Lemma 4.1, 

\G ai ,..., a Je*) - e-^ 2 / 2 | < O(N A N*0 A ) max (V^, \G au ..., am (^)\) 

= 0(N A NJ A e-' 72e2 / A ). 
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Integrating, we find 



\h\< [ G ai ,..., arn (e") - 

J\6\<t/N 

/OO / 1 \ 

e * e -* 2 e 2 /4 de = 0(N 4 N*a- 5 ) = O (-==-). 
-oo C-i »* ' 



-oo 

Further, again using Lemma 4.1, 

' ' " Jr/N ~ ~ <J(<JT/N) 2 \<tNJ 

Finally, \I 2 \ < f* /N |F ai ,..., 0m (e w )| d9. □ 

In order to verify Conjecture 1.4, it thus suffices to show that the integral 
/ T % \F ai ,..., am (^)\ d9 in Lemma 4.2 is o(^). 

Remark 4.3. For example, an estimate 

Fa 1 ,...,a m (e W ) = 0(J ¥ ), 0<9<7T, (4.1) 

is sufficient for (1.11). We conjecture that this estimate (4.1) holds when 
iV* > 6, say. Note that it does not hold for very small N*: taking 9 = ir we 
have, for even m, F ni i(— 1) = l/(^i + 1) = 1/N, and the same holds for 

Note further that even the weaker estimate 

Fa 1 ,...,a m (e W )=0(-^ ¥ ), 0<#<7T, (4.2) 

— 1/2 

would be enough to prove (1.11) with the weaker error term 0(iV* ). 

We obtain a partial proof of Conjecture 1.4 using the following lemma. 

Lemma 4.4. For a given r G (0, 1] there exists c = c(r) > such that 

\F m , n2 (e id )\ < e~ cn * (4.3) 

for m > ri2 > 1 and rj [n\ + n-i) < \9\ < it. 

More generally, for any a±, . . . ,a m and t/N < \9\ < it, 

\F ai ,..., am (e ie )\ < e~ cN *. (4.4) 
Proof. We prove first (4.3). For positive integer n define 

n 

fn(y,q) = U(i-yg i r 1 . 



3=0 

For < R < 1, we have (e.g. by Taylor expansions) e 2R < jz^, and thus 
° AR < (l-fl)^ = 1 + (l-flp • Hence, by convexity, for any real (, 

2R(i-cosC) < 1 , 2i?(l-cosC) l + ^ 2 -2ficosC ll-^l 2 

6 - + (1 - i?) 2 (1 - i?) 2 (1 - i?) 2 ' 
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and thus 

(l-ite*) -1 < (1--R) -1 exp (-12(1- cos C)). 
Consequently, by a simple trigonometric identity, for any real (f) and 9, 
f ni (Re^,e l9 )\ < (l-R)^- 1 

m s sin(ni + 1)9/2 



sin 9/2 



< 



—R\n\ + 1 — cos(0 + -^-^) ' 



The function g(#) = g n {9) := s ™™(gff , where n > 1, is an even function 
of 0; is decreasing for < 9 < n/n, as can be verified by calculating </; 
and satisfies \g(9)\ < g(ir/n) for ir/n < \9\ < ir. Further, for n > 2 and 
< \9\ < vr/n, 

9n{9) = 2 S " l( "^ ) cos(nfl/4) = 2g n/2 (9) co S (n9/4) < ncos(nfl/4) 

s "0 --«>-)• 

Let #o = T~(rii + n 2 )~ 1 < 7r/(ni + 1). For 9q < \9\ < tt we thus have 



\g ni+ i(9)\ < g ni+ i(9 ) < m + l-^; 
whence, for < R < 1, the estimate above yields 

/ m (J?e^,e i0 )| < (1 - I?)-" 1 - 1 exp (-12n^ 2 /40) . (4.5) 

Combinatorially we know that [y e q n ]f ni (y, q) is the number of partitions 
of n having at most I parts no one of which exceeds n\. As said in Re- 
mark 1.1, this equals [q n ]( n ^)F ni+ £(q). Hence, using Cauchy's integral 
formula, for any R > 0, 

dy 



3/)2 



( ni + n2 ) Fniin2 ( q ) = [y n VnAy,q) = w~* f UAv,qY 

V n l J 27T1 J\y\=R I 



'\v\=r y 

Consequently, (4.5) implies that for 9q < \9\ < tt and < R < 1, 



H2 + 1 



111 



2 ) \F nu nM\ < (1 " RY^R- 712 exp (-tfn^/40) . 

Now choose R = p := ni/{n\ + 71-2) < 1/2. By Stirling's formula, 
V+n 2 X =n(n - 1 / 2)(1 _ /9) _ m _ v _ n2 

and thus, for #0 < \9\ < it, 

\F nu nM\ < 0(nl /2 ) exp (-pn?0g/4O) = 0(n 2 /2 ) exp (-0(n 2 )) . 
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This shows (4.3) for n,2 sufficiently large. To handle the remaining finitely 
many values of ri2 we shall show: for each ri2 > 1 and r G (0, 1], there exists 
5 > such that 

\F ni ,n 2 (e w )\ < 1-8 (4.6) 
for all n\ > ri2 and t j(n\ + n.2) < \9\ < ir. To do this, we use 



l^ l5 n 2 (e ie )| = n 



:r- 



\ni + 3 



sin(m +j)6/2 



sin(j0/2) 



. (4.7) 

smj6>/2 ni+j 

Let iV = ni +ri2 and r/iV < |#| < ir. For n > n\ + 1 we have \nd\ > n\\9\ > 
N\6\/2 > r/2, and thus the estimates above show that 

\g n {0)\ <g n {r/2n) <n(l-r 2 /160) . 

Hence the final product 1I2 in (4.7) is bounded by 1 — r 2 /160 < 1. The 
product LT1 is a continuous function of 6, and equals 1 for 8 = 0; hence 
I ni| < 1 + r 2 /200 for \9\ < e, where e > is sufficiently small. (Recall that 
n,2 now is fixed.) This proves (4.6) for \9\ < e. 
For larger \6\ we use the factorization 

Fn ^ q) -{ n 2 ) (l-g)---(l-g-) • (48) 

Let < \6q\ < ir, and suppose that k > of the factors in the denominator 
of (4.8) vanish at q = q := e l9 °. Then < k < n 2 - 1, since 1 - q ^ 0. 
There are at least k factors in the numerator of (4.8) that vanish at qo (since 
F is a polynomial, and all factors have simple roots only); for q = e ld , each 
of these factors is bounded by N\q — qo\ < N\9 — 9q\ while every factor is 
bounded by 2; hence the numerator of (4.8) is O(N k \0 — 8 \ k ). Let J be 
an interval around 6q such that the denominator of (4.8) does not vanish at 
any q = e ld 7^ qo with 9 G J; then the denominator is Q(\6 — #o| fc ) f° r G J. 
Finally, the binomial coefficient in (4.8) is 0(n™ 2 ). 

Combining these estimates, we see that uniformly for 6 G J, 

K^'l - °( $v>~-w ) = 0(n " m) = 0( " rl) ' 

Since the set e < \9\ < ir may by covered by a finite number of such interval 
\Fn!,n 2 ( el0 )\ = 0(n ] ~ 1 ) uniformly for e < \6\ < ir. Consequently (4.6) 
holds for all such 9 if n\ is sufficiently large. 

It remains to verify (4.6) for each fixed 712 and a finite number of n\\ in 
other words, that for each ri2 > 1 and ri\ > 112, there exists 5 > such that 
(4.6) holds. To see this, note that the events M nitU2 = and M nitU2 = 1 
both have positive probability. It follows that |-f 1 ni,n 2 ( el6l )l < 1 f° r every 6 
with < |#| < 7T, and (4.6) follows. This completes the proof of (4.3). 
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To prove (4.4), we assume as we may that a\ > ■ ■ ■ > a m and use the 

factorization (1.8). Let J be the first index such that 02 H + aj > iV*/2. 

For j > J, then A,_i + aj = Aj > Aj > ai + N*/2 > N/2, and thus 
Aj\0\ > N\6\/2 > r/2; hence (4.3) yields 

\F Aj _ u aM 6 )\< e ~ C{T/2)ai - 

We thus obtain from (1.8), since each F nitn2 is a probability generating 
function and thus is bounded by 1 on the unit circle, 

m m 

\F ai ,...,a m (e*)\ = J] \F Aj . u aM 6 )\ < II e- C(T/2)a > < e-^ N *'\ 

3=2 j=J 

because *Y^j=j a i — ^*/2- This proves (4.4) (redefining c(r)). □ 

Theorem 4.5. There exists a positive constant c such that for every C , the 
following is true. Uniformly for all a\, . . . , a rn such that a* < Ce cN * and all 
integers k, 

"- = fc > = ^^"-" ,1/,2 ' 2,+0 (^))^ < 4 - 9 > 

Proof. Let c\ = c(r) be the constant in Lemma 4.4. Then, Lemmas 4.2 
and 4.4 yield 

- = *) = ^{' Hk '^' ) + °(w, + ' e ~ nN -)y 

For any fixed c < c\ we have, using Lemma 3.1, aN*e~ ClN * = 0(Ne~ cN *) 
and thus 

The result follows, since Ne~ cN * = a*e~ cN * +N*e~ cN * = a*e~ cN * +0(1). □ 

4.1. Log-concavity. Let us review the proof of Theorem 4.5 with the in- 
tention of greater accuracy. The goal is to prove log-concavity in some range. 
For concreteness, let a = (n,n). Then a 2 is of order n 3 , and for sufficient 
accuracy we take the Taylor series in the exponent of (3.2) out to 0(9 10 ). 
This yields, for some polynomials pk(n) of degree k + 1, 

F ntn (e ie ) = exp(i^ - a 2 6 2 /2 + p 4 (n)9 4 + p 6 (n)9 6 + p 8 (n)9 8 + 0(n n 9 10 )) 
= e ^- 2e2 /2 (! + pA{n)9 4 + pe{n)e 6 + pg{n)e 8 

+ \pl{n)6* +p A {n)p Q (n)e w + \pl(n)9 12 + 0(n n 6 10 )) 

Arguing as in the proof of Lemma 4.2 but using this estimate instead of 
(3.5) for \9\ < t/N, one easily obtains, after the substitution 8 = t/a, for 
any k and with x := (k — n)/o, 

n»(M„ = *j = -L £ e - 2 ' 2 -"-(i + ^ + ■ ■ ■ + fS* 12 ) I + om-V). 
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Letting <p(x) := (27r) _1 / 2 e _a;2 / 2 denote the normal density function, and <pW 
its derivatives, we obtain by Fourier inversion 

P(M B>B = k) = a-' ( V (x) + ^^\x) + ■■■ + f^ (12) (*) + 0(n~ 4 )) 

1 e -xV2(! + QM + (n- 4 )), (4.10) 
lira 



where Q(n, x) is a' 12 times a certain polynomial in n and x of degree 17 in 
n; thus for cc = 0(1) we have Q(n, x) = 0(n _1 ) and similarly, for derivatives 
with respect to x, Q'(n,x) = 0(n _1 ) and Q"(n,x) = 0(n~ l ). (Q(n,x) can 
easily be computed explicitly using computer algebra, but we do not have 
to do it.) 

Replacing k by k ± 1 in (4.10) we find, for x = O(l), 

P(M„, n = k ± 1) = -J^ e -(*±0 2 / 2 (l + g (n> x ) ± CT -iQ'( n) a) + 0(n" 4 )) , 
V 2tto~ 



and thus 

P(M n ,„ = A:-l)P(M nin = A: + l 
1 



2 ^ 2 e- x {(l + Q(n,x)) 2 -a~ 2 Q'(n,x) 2 + 0(n- 4 )). 



e~^ 2 ¥(M nyn = k) 2 (l + 0(n- 4 )). 
Hence, for x = O(l), i.e., k = n + 0(a), 
P(M n , n = fc) 2 - P(M n , n = k - 1) P(M n> „ = fc + 1) 

= (a- 2 + 0(n~ 4 )) P(M n , n = fc) 2 = ^e- 2 (1 + 0(n- x )) . 

(4.11) 

In particular, this is positive for large n. This gives: 

Theorem 4.6 (A log-concavity result). For each constant C we have no 
such that for n > uq and \ j — (i\ < Ca 



Cj > Cj-lCj +1 , 



where 



q 

We note that the "mysterious" numbers appearing in our earlier table for 
the choice j = n 2 /2 — 1 are asymptotically 



1 {2n\ 2 18 {2n\ 2 18 



27TO" 4 \ n 7rn 6 \ n it 2 



n~'2 



7Mn 



Remark 4.7. This argument for log-concavity in the central region does 
not use any special properties of the distribution; although we needed sev- 
eral terms in the asymptotic expansion above, it was only to see that they 
are sufficiently smooth, and the main term in the final result (4.11) comes 
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from the main term e~ x ' ''/ 2 / '(y/2n<r) in (4.9). What we have shown is just 
that the convergence to the log-concave Gaussian function in the local limit 
theorem is sufficiently regular for the log-concavity of the limit to transfer 
to F(M n>n = k) for k = fi + 0(a) and sufficiently large n. 

5. Final comments 

Suppose that iV* -/-> oo. We may, as usual, assume that a\ > ■ ■ ■ > a m . 
By considering a subsequence (if necessary), we may assume that := 
N — a* = a>2 + • • • + a m is a constant; this entails that m is bounded, so by 
again considering a subsequence, we may assume that m and ci2, . . . , a m are 
constant. We thus study the case when a\ — ► 00 with fixed 02, . . . , a m . 

In this case, the number of inversions between indices 2, ... ,m is 0(1), 
which is asymptotically negligible. Ignoring these, we can thus consider the 
random word as iV* letters 2, . . . ,m inserted in a± l's, and the number of 
inversions is the sum of their positions, counted from the end. It follows 
easily, either probabilistically or by calculating the characteristic function 
from (1.6), that M aij .„ >am /N, or equivalently M aij ... )0m /ai, converges in dis- 
tribution to the sum Uj of independent random variables Uj with 
the uniform distribution on [0,1]. Equivalently, since a 2 ~ n\N*/l2 ~ 

iv 2 ivyi2, 




where — ► denotes convergence in distribution. This limit is clearly not 
normal for any finite N*. (However, its distribution is close to standard 
normal for large N*. Note that it is normalized to mean and variance 1.) 
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